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SUMMARY 

This  report  discusses  two  methods  of  liquid  phase 
epitaxial  (LPE)  growth  applied  to  the  compound  semi¬ 
conductor  gallium  arsenide,  and  develops  a  mathematical 
description  for  each  using  a  diffusion-limited  model. 

The  properties  of  the  epitaxial  layers  grown  under 
various  experimental  conditions  are  discussed,  and  the 
observed  growth  rates  compared  with  those  predicted  by  the 
mathematical  models.  From  the  differences  observed,  it  is 
concluded  that  processes  additional  to  simple  diffusion  of 
solute  through  the  melt  are  occurring.  Convective  mixing 
caused  by  horizontal  temperature  gradients  within  the  melt 
is  suggested  as  an  explanation  for  these  differences. 
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LIST  OP  FIGURES 

1.  Schematic  of  apparatus  used  for  liquid  phase  epitaxial  growth  of  a  single 
layer. 

(a)  Sliding  boat  used  for  the  Temperature  Transient  Technique. 

(b)  Crucible  arrangement  used  for  the  Steady  State  Method. 

2.  (a)  Various  surface  morphologies  observed  on  Sample  8. 

(Magnification  140X) 

(Scale:  71  ju  /cm  ) 

(b)  Various  surface  morphologies  observed  on  Sample  9. 

(Magnification  140X) 

(Scale:  71  /i/cm) 

3.  Irregular  "mesa"  structures  formed  when  the  saturated  melt  does  not  wet 
the  substrate  surface  uniformly. 

(a)  Surface  features  (Scale:  28.6  /i/cm) 

(b)  Cross-sectional  view  of  cleaved  sample. 

(Scale:  13.8  /i/cm) 

4.  Highly  uniform  junction  formed  on  a  GaAs  substrate. 

(Scale:  11.4  /i/cm) 

5.  Irregular  junction  formed  after  attempts  to  dope  the  epitaxial  layer 
with  silicon.  (Scale:  66  /i/cm) 

6.  Typical  example  of  the  variation  in  thickness  of  the  epitaxial  layers 
grown  using  the  Steady  State  Method.  The  sample  shown  here  was  grown 
over  a  period  of  two  hours,  at  a  nominal  temperature  of  800°C  and  with 
a  temperature  differential  of  50  C  applied  across  the  furnace. 


1 


ERL-0078-TR 


1 .  INTRODUCTION 

Epitaxial  growth  techniques  provide  excellent  control  over  the  purity 
and  the  crystalline  perfection  of  many  semiconducting  materials.  With 
GaAs  for  example  the  semiconductor  injection  laser  is  usually  made  from 
epitaxially  grown  material,  because  the  very  low  defect  density  and  the 
background  impurity  level  required  to  achieve  efficient  radiative 
recombination  are  difficult  to  obtain  by  non-epitaxial  methods.  The 
phenomenon  of  epitaxial  growth  is  of  considerable  scientific  interest  in 
its  own  right,  but  in  addition  the  semiconductor  lasers  produced  by 
epitaxy  are  small  and  efficient  and  of  increasing  importance  in  the  fields 
of  communications  and  surveillance.  It  is  these  considerations  which 
have  led  to  the  work  on  epitaxial  growth  of  GaAs  to  be  described  in  this 
technical  report. 

Chemical  vapour  phase  epitaxy  and  liquid  phase  epitaxy  are  both  used 
extensively  for  the  growth  of  semiconductors.  However  for  the  III-V 
compounds  gas  phase  epitaxy  tends  to  give  poor  compositional  control  in 
the  grown  layers  because  of  difficulties  in  accurately  controlling  the 
gas  phase  compositions  and  because  the  group  V  elements  dissolve  in  the 
group  III  source  material.  Thus  for  these  compounds  liquid  phase  epitaxy 
is  preferred. 

In  this  technical  report  we  will  examine  some  theoretical  and  practical 
aspects  of  two  methods  of  liquid  phase  epitaxial  growth,  the  temperature 
transient  technique  and  the  steady  state  differential  technique.  Because 
Laser  Group  possessed  no  facilities  for  epitaxial  growth  prior  to  1976,  the 
experiments  described  in  this  report  were  carried  out  at  the  Australian  Post 
Office  (now  Telecom  Australia)  Research  Laboratories  in  Melbourne  from  April 
to  July  1974,  and  in  the  Electrical  Engineering  Department  of  Adelaide 
University  from  November  1974  to  July  1975.  At  both  centres  the  work  was 
necessarily  restricted  by  equipment  which  had  been  designed  for  other  purposes. 
However,  the  experience  gained  during  these  periods,  has  allowed  the  design 
of  a  system  for  liquid  phase  epitaxy  which  has  now  been  set  up  at  DRCS. 


2.  METHODS  FOR  LIQUID  PHASE  EPITAXIAL  GROWTH 

Because  the  solubility  of  a  dilute  constituent  such  as  arsenic  in  a 
liquid  solvent  such  as  gallium  decreases  with  decreasing  temperature,  then  if 
the  temperature  of  an  initially  saturated  solution  is  reduced, 
recrystallisation  of  gallium  arsenide  will  take  place.  In  the  liquid  phase 
epitaxial  process,  conditions  are  controlled  so  that  the  recrystallised 
material  forms  as  a  uniform  monocrystalline  deposit  on  a  single-crystal 
substrate  wafer. 

In  the  temperature  transient  technique , (TTT) ,  used  at  the  Telecom  Australia 
Research  Laboratories,  molten  gallium  doped  with  selected  impurities  is 
saturated  with  arsenic  using  an  undoped  gallium  arsenide  source,  and  then 
placed  in  contact  with  a  gallium  arsenide  single  ciystal  substrate.  When 
this  solution  is  cooled  the  melt  becomes  supersaturated,  which  causes 
precipitation  of  doped  gallium  arsenide  and  its  subsequent  nucleation  in 
single  crystal  form  on  the  gallium  arsenide  substrate.  The  rate  of  growth 
of  the  epitaxial  layer  depends  upon  the  initial  temperature  and  the  rate  of 
temperature  decrease.  Growth  is  terminated  after  the  desired  time  interval 
by  mechanically  pushing  the  substrate  out  of  contact  with  the  melt.  A  slight 
variation  of  the  technique  used  here  allows  the  growth  of  multiple  layers 
by  pushing  the  substrate  successively  into  contact  with  different  melts, 
but  for  the  work  in  this  report  only  single  layers  were  grown. 
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Because  this  type  of  epitaxial  growth  takes  place  at  temperatures  well 
below  the  melting  point  of  the  III-V  compounds,  high  residual  purity  is 
generally  obtained  since  the  major  impurities  have  low  distribution 
coefficients  at  these  temperatures.  However,  the  effective  impurity 
concentration  changes  with  both  temperature  and  growth  rate,  resulting  in  a 
spatially  varying  dopant  concentration  as  the  solution  is  cooled.  For  very 
thin  layers  this  effect  may  be  unimportant;  however,  for  thick  layers  the 
variation  in  dopant  concentration  across  the  layer  may  be  unacceptable. 

In  this  situation  it  may  be  advantageous  to  use,  as  was  done  at  Adelaide 
University,  the  steady  state  differential  method  (SSM)  where  growth  takes 
place  at  a  constant  temperature.  A  temperature  gradient  established  within 
the  melt  provides  the  driving  force  for  diffusion  of  arsenic  from  the  source 
gallium  arsenide  through  the  melt,  to  be  deposited  on  the  cooler  substrate 
as  an  epitaxial  layer  of  gallium  arsenide.  Although  thick  layers  with  uniform 
dopant  concentration  can  be  grown  using  SSM,  the  thickness  of  these  layers 
cannot  be  as  accurately  controlled  as  in  the  TTT.  This  lack  of  precision 
is  partly  due  to  the  greater  number  of  variables  which  must  be  accurately 
controlled  in  the  SSM,  but  in  addition,  the  actual  growth  processes  in  the 
SSM  are  less  well  understood  than  in  growth  by  TTT.  This  fact  will  become 
evident  from  the  work  discussed  in  this  report. 

In  either  method  of  LPE,  slow  growth  rates  encourage  the  homogenous 
incorporation  of  dopant  impurities  including  those  with  distribution 
coefficients  greatly  different  from  unity.  With  sufficiently  slow  growth 
rates  the  melt  composition  is  essentially  constant  right  up  to  the  growing 
film-melt  interface,  whereas,  for  high  growth  rates,  growth  may  occur  from 
a  melt  whose  local  composition  can  be  very  different  from  its  average  value. 

In  addition  the  temperatures  used  for  LPE  growth  are  several  hundred  degrees 
below  the  melting  point  of  the  solid  material  being  formed,  so  that  diffusion 
from  the  substrate  into  the  growing  layer  is  minimised.  This  reduced 
temperature  (compared  to  the  usual  bulk  crystal  growth  methods  such  as 
Czochralski)  decreases  the  level  of  contamination  introduced  into  the  melt 
from  the  growth  system,  and  results  in  lower  background  impurity  levels  in 
the  growing  layer. 


3.  MATHEMATICAL  FORMULATION  OF  LPE  GROWTH 

By  making  several  simplifying  assumptions  about  the  growth  processes,  it 
is  possible  to  derive  mathematical  models  for  the  height  of  epitaxially 
grown  layers.  The  validity  of  the  a.ssumptions  can  then  be  tested  by  comparing 
the  predicted  layer  heights  with  those  observed  experimentally.  The 
assumptions  used  in  deriving  the  present  model  are  listed  below: 

(i)  Transport  of  arsenic  to  the  substrate  takes  place  by  diffusion 
through  the  melt  from  a  source  having  uniform  composition,  to  a 
sink  (the  growth  surface)  whose  concentration  at  any  time  is 
determined  by  the  temperature  and  therefore  depends  upon  the 
particular  growth  technique  employed. 

(ii)  After  the  substrate  is  placed  in  the  solution,  the  liquid  melt 
and  solid  substrate  are  considered  to  be  in  equilibrium,  with  the 
concentration  of  arsenic  at  the  growth  interface  determined  from 
the  liquidus  properties  appropriate  to  the  temperature  of  the 
melt. 

(iii)  Excess  solute  arsenic  within  several  diffusion  lengths  of  the 
growing  interface  is  removed  from  solution  only  by  epitaxial 
deposition  on  the  substrate,  and  the  presence  of  the  substrate 
causes  no  edge  effect  perturbations. 
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(iv)  No  reerystallisation  of  gallium  arsenide  occurs  in  the  melt. 

This  assumption  is  realistic  since  it  is  known  that  the  dilute 
solutions  used  for  LPE  growth  of  gallium  arsenide  will  withstand 
supercooling  far  in  excess  of  that  used  for  LPE,  without 
spontaneous  nucleation  occurring. 

(v)  There  is  no  crystallisation  of  arsenic  on  the  melt  surface,  and 
no  loss  of  arsenic  through  evaporation  into  the  hydrogen 
atmosphere  which  is  the  growth  environment. 

(vi)  The  diffusion  coefficient  and  the  slope  of  the  liquidus 
(solubility)  curve  may  be  considered  constant  over  the  small 
temperature  interval  of  each  growth  run. 

(vii)  Temperature  uniformity  exists  within  the  solution  in  that  there  are 
no  unplanned  thermal  gradients  which  could  cause  kinetic 
processes  (such  as  convection)  and  perturb  the  uniformity  of 
solution  at  the  growth  interface., 


Very  little  fundamental  data  is  available  regarding  this  last  assumption, 
but  as  will  become  evident  later  in  this  report,  such  factors  may  have  a 
very  great  influence  on  the  growth  rates. 

The  temperature  transient  growth  technique  is  the  experimental  method  most 
frequently  applied  in  overseas  device  technology,  and  some  quite  elaborate 
mathematical  models  have  been  developed  (references  5  to  1 1 ) .  In  the  present 
work  we  have  used  a  simplified  approach  based  directly  on  the  variation  of 
arsenic  solubility  with  temperature.  Such  an  approach  was  suggested  by 
Teede  (reference  12),  but  his  formulation  has  been  somewhat  modified  by 
the  present  author.  After  developing  a  model  for  temperature  transient  growth, 
we  will  extend  this  approach  to  develop  a  model  for  the  steady  state 
temperature  gradient  technique.  This  latter  model  was  derived  from 
discussions  with  McRae  (reference  13). 

3.1  Temperature  Transient  Technique  (TTT) 

For  simplicity  in  the  following  formulations,  the  chemical  symbols 
Ga,  As  and  GaAs  will  be  used  for  gallium,  arsenic  and  gallium  arsenide 
respectively. 

With  the  notation  that 


XAs(T> 


solubility  (atom  fraction)  of  As  in  Ga  at  temperature  T 


nGa(melt) 


number  of  atoms  of  Ga  solution  added  to  the  source  from 
the  melt 


n~  ,  (source, 
GaAs 


number  of  atoms  of  GaAs  dissolved  from  the  source 
material 

nG&( source)  +  n^g( source) 


nu  (melt) 
Ga 


mass  of  Ga  solution  added 


mGaAs(T) 


M 


Ga 


M, 


GaAs 


=  mass  of  source  GaAs  in  solution  at  temperature  T 
=  molecular  weight  of  Ga  (=69.72) 

=  molecular  weight  of  GaAs  (=144.63) 
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then  (t)  =  solubility  (atoEi  fraction)  of  As  in  Ga 


n^R (source) 
n,  (source)  + 

Aft 


n~  (source)  + 
Ga 


n„  (melt) 
Ga 


Now  n.  (source) 

A  Q 


therefore 


nn  .  (source)  =  n~  (melt) 
GaAs  Ga 


V(T) 


1-2Xls(T) 


The  mass  of  GaAs  dissolved  from  the  source  material  and  in  solution 
at  temperature  T^  is  then 


m„  .  (T  ) 
GaAs  o 


M 


GaAs .  mr  (melt) 


V(T) 


1-2XAs(T) 


(1) 


The  amount  of  GaAs  deposited  from  the  source  onto  the  substrate  during 
a  temperature  reduction  from  T^  to  T^  ,  GaAs^ep^’ 


mGaAs^Dep^  =  MGaAs  . 

MGa 


XA  (T  )  X,  (T. ) 

— A§ — Q — ..  ~  - — 1 — - 


Using  the  abbreviation  X^ 


(T)  =  XAa^ 

’-2VT) 


then 


mGaAs(Dep)  =  2.074  mGa(melt) 


The  density  of  GaAs  is  5.32  g/cc,  so  that  the  height,  h,  of  the 
epitaxial  layer  grown  on  a  GaAs  substrate  of  area  A  is 


h  =  •  mGa(melt) 


x;  (t  )  -  x!  (t. ) 

As  o  As  1 


(2) 


Layers  grown  at  the  Telecom  Australia  Research  Laboratories  using  the 

TTT  on  substrate  areas  of  approximately  0.6cm  and  solvent  masses  mGa(melt) 

of  approximately  1  gram,  over  temperature  intervals  (T  -T  )  of  50°C  and 

100°C  and  initial  temperatures  Tq  of  700°C  and  800°C,  were  observed  to  grow 

to  approximately  60$  of  the  magnitudes  predicted  by  equation  (2)  when  the 
liquidus  data  of  Hall  (reference  14)  was  used.  This  agreement  is  regarded 
as  satisfactory  in  view  of  the  greatly  simplified  approach  taken  in  deriving 
equation  2.  Several  quite  complex  theoretical  formulations  for  temperature 
transient  growth  have  shown  no  better  agreement  with  experimental  observations 
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3.2  Steady  State  Method 

We  will  now  consider  a  model  suitable  for  use  with  the  Steady  State 
Differential  Method,  in  which  a  temperature  differential  AT  imposed 
across  a  solution  causes  a  corresponding  equilibrium  solute  concentration 
differential  AC  between  source  and  substrate.  In  the  absence  of  information 
on  the  fluid  dynamics  of  the  melt  solution,  we  will  formulate  an  equation 
under  the  assumption  that  solute  diffusion  is  the  process  which  limits  the 
rate  of  growth  of  the  epilayer. 

In  addition  to  the  notation  expressed  in  3.1,  we  will  use  the  following 
in  this  section 

=  solute  (source  GaAs)  concentration  in  the  Ga  melt,  expressed  in 
terms  of  the  mass  of  dissolved  GaAs  per  unit  volume  of  Ga  melt. 

=  diffusion  coefficient  of  As  under  a  concentration  gradient  in 
the  Ga  melt. 

=  temperature  gradient  applied  across  the  Ga  melt. 

=  solute  concentration  gradient  through  the  melt  caused  by  the 
applied  temperature  gradient. 

=  variation  in  solute  concentration  with  temperature. 

=  density  of  GaAs  =  5.32  g/cc. 

=  density  of  Ga  =  5.95  g/cc. 

From  the  law  of  diffusion  we  obtain  the  rate  of  mass  deposition  per  unit 
area  A  to  be 


C 

D 

dT 

dx 

dC 

7GT 

dC 

dT 

pGaAs 

pGa 


1  £_ 

A  dt 


m, 


GaAs 


dll 

so  that  the  rate  of  growth,  of  the  deposited  layer  becomes 


dh 

dt 


-D 

pGaAs 


dC 

dx 


D _  dC_  dT 

PGaAs  '  dT  ‘  dz 


Now  from  equation  (l  ) 


in  3.1 ,  we  know  that 


mGaAs 


(T) 


GaAs.  . 
Ga 


m^a(melt) 


(T)  . 
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Using  this  expression,  and  the  definition  of  concentration  in  terms  of 
dissolved  mass  per  unit  volume  of  melt,  the  concentration  C  of  the 
solute  near  the  source  is  then 


C 


GaAs 

Ga 


(T  ) 

source 


and  the  concentration  difference  AC 
Tsource  and  substrate  at  temperature 


between  source 

T  ,  is  then 
sub 


at  temperature 


AC 


M 


GaAs 


M, 


Ga 


Ga 


x'  (T 

As  source 


)  -X1  (T  ) 

As  sub ' 


=  M. 


GaAs 


M 


Ga 


o  AX1  (T  -*  T  ) 

As  source  sub 


(3) 


Since  the  slope  of  the  liquidus  curve  is  constant  over  a  small 
temperature  interval  AT,  then 


dC  Ac  p  MGaAs  (  AXAS  ^ 

dT  AT  Ga  H  V  AT  / 

Ga 

If  we  can  also  assume  that  the  temperature  gradient  within  the  solution. 
AT 

_ §.  >  1S  constant  over  the  distance  between  source  and  substrate,  then 

Ax 


AT 


s  =  dT 


A  x 


dx 


so  that 


dh  _  Ga 


M 


GaAs 


dt 


VL 


GaAs  Ga 


h) 

at 


sub 


(4) 


When  the  diffusion  coefficient  D  is  expressed  in  terms  of  cm^/sec,  and 
tne  solution  temperature  gradient  _ s  in  C/cm,  then  the  expression 


for  rate  of  growth  of  the  epilayer  dh  becomes 

dt 


at  -  2-32D 


"at 

Ax 


A 

slope  of  curve  X'  (t)  -  T 

AS 


(cm/sec)  or,  of 


IT 


sub 


more  practical  use. 


£jT  =  8.35  x  1 07  D 


at, 

Ax' 


slope  of  curve  X,  (t)  -  T 
As 


sub 


{u m/hour) 
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Long  et  al  (reference  15)  have  developed,  a  more  complicated  diffusion 
equation: 

Ac  (substrate-solution  interface) 

AC  (source-solution  interface) 

which  under  similar  ass\imptions  to  equation  (4)  simplifies  to 


dh  _  /  D_\  log 
dt  \  Ax/  ' 


dh 

dt 


D 


sub 


a 


exp 


y 

(-  +  V 

T  v 
sub 


(5) 


where  C  is  the  concentration  of  As  in  solid  GaAs  and  a  ,  /3  ,  y  are 

S  Arp  \ 

constants.  When  the  same  values  for  D  and  1 _ s  are  inserted  in  both 

\  Ax/  dh 

equations  (4)  and  (5),  the  calculated  values  for  dt  are  substantially 
the  same.  Over  the  range  of  interest,  the  difference  is  about  5$. 

To  calculate  values  for  ^  with  which  to  compare  experimental 

observations  we  must  first  determine  values  for  the  diffusion 

AT 

coefficient  D,  the  temperature  gradient  _ s  ,  and  the  slope  of  the 

Ax 

liquidus  curve  at  each  temperature  of  interest.  Such  values  are 
discussed  in  Appendix  1 ,  and  the  values  shown  in  Table  1  are  considered 
appropriate  to  the  present  situation. 


TABLE  1 . 


T  (°C) 

D  (cm2/ sec) 

Slope 

4 (T)  -  T_ 

T  V 

sub 

at  0  , 

_ s  C/cm 

A  x 

700 

800 

3  x  10-5 

5  x  10"5 

-4 

0.90  x  10 

2.78  x  10~4 

4.2 

4.6 

Using  these  values  with  equation(4),  grown  layer  heights  of  1  .Op /hour 
at  700°C  and  5.3  M/hour  at  800°C  are  predicted  for  the  SSM.  The  rates 
observed,  about  5  p/hour  and  15  p/hour  at  700  C  and  800  C  respectively, 
are  much  larger  than  the  values  predicted  on  the  basis  of  diffusion 
limited  growth.  However,  as  reported  in  reference  16,  convective  forces 
caused  by  temperature  non— uniformities  in  the  furnace  were  found  to  have 
contributed  to  the  epitaxial  growth  and  greatly  increased  the  deposition 
rate  of  GaAs  from  the  saturated  melt. 


4.  EXPERIMENTAL  DETAILS 

The  various  techniques  for  liquid  phase  epitaxial  growth  allow  the 
formation  of  highly  crystalline  material  with  low  defect  and  background 
impurity  concentrations,  together  with  the  ability  to  introduce  donor  and 
acceptor  impurities  in  a  strictly  controlled  manner.  To  achieve  these 
properties,  the  individual  aspects  of  the  growth  process  must  be  carefully 
considered  and  optimised.  The  major  areas  requiring  detailed  consideration 


are: 
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(i) 

temperature  control 

(ii) 

growth  vessel 

(iii) 

growth  environment 

(iv) 

substrate  material  preparation 

(v) 

source  material  preparation 

(vi) 

melt  preparation 

(vii) 

impurity  dopants 

(viii) 

loading  procedure. 

4.1  Temperature .control  . 

In  both  growth  methods,  stability  and  accurate  control  of  the  melt 
temperatures  are  imperative.  The  SSM  used  a  split  furnace  whose  upper 
element  was  maintained  at  a  lower  temperature  than  the  bottom  element, 
creating  a  temperature  gradient  across  the  furnace  and  thus  across  the 
melt  solution.  The  temperature  gradient  and  the  actual  temperature  in  the 
melt  should  be  constant  over  at  least  the  region  where  the  melt  contacts 
the  substrate.  The  great  importance  of  this  requirement  is  demonstrated 
in  later  sections.  The  furnace  used  in  the  TTT  had  a  uniform  temperature 
throughout  the  operating  volume,  and  this  temperature  was  accurately 
controlled  to  produce  a  steady  cooling  rate  in  the  furnace  over  a  period 
of  several  hours.  In  the  experiments  reported  here,  furnace  temperatures 
in  the  range  of  600°  -  900°C  were  used,  with  applied  temperature 
differentials  of  0  ,50  or  100°C,  and  various  cooling  rates  in  the  range 
of  15°/hour  to  50°C  hour.  The  various  furnace  temperatures  were  set  using 
commercial  controllers  whifch  were  supplied  with  the  respective  furnaces, 
while  the  cooling  rates  in  the  TTT  were  controlled  by  an  electronic  ramp 
generator.  Temperatures  were  monitored  in  all  cases  using  chrome 1-alumel 
thermocouples,  although  it  has  since  been  realised  that  this  was  not  the 
optimum  choice. 

4.2  Growth  vessel 

A  cylindrical  quartz  tube  within  the  furnace  contained  the  growth 
vessel,  whose  background  contamination  had  previously  been  reduced  to  a 
low  level  by  vacuum  baking  at  high  temperatures  followed  by  baking  for  an 
extended  period  in  a  flowing  hydrogen  atmosphere  while  loaded  with 
gallium.  For  the  TTT  a  machined  pyrolytic  graphite  sliding  boat  with 
several  wells  or  bins  containing  the  respective  melts  was  used  for  multiple 
layer  growths.  The  substrate  wafer  was  moved  from  bin  to  bin  at  certain 
times  in  the  temperature  cycle,  and  at  the  completion  of  growth  was  moved 
out  of  contact  with  any  melt.  The  simpler  boat  which  was  used  for  the 
growth  of  the  single  layers  reported  here  is  illustrated  schematically 
figure  1  (a). 

The  vessel  in  use  at  Adelaide  University  for  the  SSM  is  illustrated 
in  figure  l(b).  This  consisted  of  a  vitreous  carbon  crucible  in  which 
the  source  block  was  held  down  by  a  graphite  peg.  The  substrate,  typically 
5  mm  x  10  mm,  was  held  by  suction  onto  a  pyrolytic  graphite  holder. 

Movement  of  the  substrate  and  holder  was  manually  controlled  from  outside 
the  furnace  by  a  long  quartz  rod,  and  growth  was  initiated  (or  terminated) 
by  rotating  the  substrate  holder  down  into  (or  up  out  of)  the  melt 
solution. 

4.3  Growth  environment 

The  crucible  environment  must  be  controlled  during  the  growth  process 
to  avoid  oxygen  contamination. 
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Any  oxygen  introduced  during  the  loading  operations  or  from  small  air 
leaks  in  the  system  influences  the  properties  of  the  growing  layer 
(references  17,  18,  19).  Adsorbed  oxygen  on  the  surface  of  the  carbon 
growth  vessel  or  the  walls  of  the  quartz  furnace  tube  is  driven  off  at 
the  temperatures  used  during  the  growth  cycle  and  rapidly  forms  a  layer 
of  Ga~0„  on  the  melt  surface  (reference  20 ).  Such  a  layer  inhibits 

'wetting'  of  the  GaAs  substrate  by  the  liquid  Ga  melt,  and  causes  areas 
of  uneven  growth  upon  the  substrate.  It  is  therefore  advantageous  for 
the  growth  to  be  carried  out  in  a  reducing  atmosphere,  such  as  hydrogen, 
which  has  been  previously  been  purified  by  passage  through  a  heated 
palladium  element  and  a  liquid  nitrogen  trap. 

At  a  melt  temgerature  of  the  order  of  800  G  any  surface  trioxide, 

Ga203  (B.P.  2500  C),  will  react  with  hydrogen  to  form  the  suboxide,  Ga^O 

(B.P.  700°C),  (reference  20)  which  will  then  be  vaporised  and  removed 
by  the  hydrogen  steam  (which  is  flowing  at  a  fixed  rate  within  the  range 
of  100-200  cc/minute) 

Ga203(s)  +  2H2  (g)  - »  Ga20  (g)  +  2^0  (g) 

However,  a  trioxide  layer  may  also  be  formed  between  the  gallium  melt 
and  the  carbon  walls  of  the  growth  vessel.  At  high  temperatures  the 
trioxide  is  converted  to  a  gaseous  suboxide  by  reacting  with  the  gallium 
melt. 

Ga203  (s)  +  4Ga  (0  — *3  Ga20  (g) 

Complete  removal  of  this  Ga20  then  depends  on  the  suboxide  dissolving  in 

the  melt  and  vaporising  at  the  melt-hydrogen  surface.  This  will  be  a 
slow  process  because  of  its  low  solubility  in  the  melt.  Extended 
periods  for  equilibrium  between  the  melt,  the  source  material  and  the 
growth  vessel  are  thus  necessary  before  the  substrate  is  brought  into 
position  and  the  growth  cycle  commenced.  The  long  equilibrium  time  also 
minimises  impurity  levels,  since  impurities  from  the  source  are  largely 
removed  by  extraction  into  the  gallium  melt  and  vaporisation  into  the 
hydrogen  stream. 

The  volatile  oxides  and  impurities  driven  from  the  melt  could  react 
with  the  substrate  surface  and  cause  an  irregular  masking  effect.  This 
would  later  result  in  uneven  contact  of  the  gallium  melt  on  the  substrate 
surface,  and  cause  consequent  non-uniform  epilayer  growth.  It  was  found 
that  more  uniform  junctions  and  better  epilayers  were  formed  if  contact 
of  the  substrate  with  the  volatile  impurities  was  avoided,  and  this  could 
be  achieved  in  both  growth  systems  by  suitable  positioning  of  the 
substrate  for  the  long  period  before  it  was  brought  into  the  growth 
position. 

4.4  Substrate  material  preparation 

Cutting  wafers  from  a  bulk  GaAs  ingot  with  a  diamond  saw  may  produce 
work  damage  extending  into  the  wafer  at  least  25  /im  from  each  sawed 
surface  (reference  21 ).  Layers  grown  directly  on  such  damaged  surfaces 
exhibit  gross  surface  and  structural  irregularities  and  so  the  damage 
was  removed  from  the  substrate  by  mechanical  and  chemical  polishing  before 
epitaxial  growth  was  commenced. 

For  these  experiments,  each  wafer  was  then  cut  into  individual 
substrates  5  mm  x  10  mm  using  a  wire  saw.  Three  or  four  substrates  were 
positioned  equidistant  around  the  perimeter  of  a  flat  circular  glass 
plate  2-g-"  diameter  and  -g-"  thick  using  resin  wax  (glycol  pthalate)  as 
adhesive. 
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The  substrates  on  the  plate  were  mechanically  lapped  using  3200  mesh 
carborundum  suspended  in  glycerine,  followed  by  careful  washing  in 
distilled  water.  This  disc  was  then  placed  under  a  similar  disc 
la  block  y"  thick  was  found  necessary  to  produce  sufficient  pressure, 

2 

30-60  g/ cm  ,  for  flat  and  regular  polishing),  and  chemical  polishing 
undertaken  by  mechanical  motion  on  a  Pellon  Microcloth.  A  1  :40 
solution  of  "White  King"  commercial  bleach  in  distilled  water  was  used 
as  the  source  of  hypochlorite  for  the  polishing  liquid.  Polishing  action 
occurred  through  the  formation  on  the  substrate  surface  of  gallium  oxide, 
which  was  then  removed  by  the  Pellon  pad.  No  abrasive  was  used  in  the 
chemical  polishing  technique,  so  that  work  damage  was  held  to  a  minimum. 
For  this  reason,  unless  saw  damage  on  the  original  wafer  was  very  bad, 
it  was  preferable  to  omit  the  mechanical  polishing  stage  and  proceed 
directly  to  chemical  polishing.  Because  dust  particles  produce  fine 
scratches  on  the  substrate  surface,  it  was  important  for  all  polishing 
to  be  carried  out  in  a  clean  environment.  After  polishing,  the  substrates 
were  removed  from  the  glass  disc  by  heating  the  resin  wax,  and  cleaned 
with  acetone  and  alcohol  before  storage  under  methanol  in  individual 
containers. 

Immediately  prior  to  loading  into  the  growth  system,  the  substrate  was 
ultrasonically  cleaned  sequentially  in  acetone  and  alcohol,  etched  in  a 
3:1  si  mixture  of  H^SO^sH^O^ sH^O  for  5-10  seconds  to  remove  the  thin 

surface  oxide  layer  and  thus  expose  a  fresh  clean  surface  of  GaAs  for  the 
growth  interface,  then  rinsed  in  distilled  deionized  water  and 
ultrasonically  cleaned  in  methanol  or  isopropyl  alcohol.  The  figures 
quoted  above  for  etchant  composition  and  for  the  time  of  immersion  were 
decided  upon  as  optimum,  as  it  was  found  that  much  longer  in  the  etch 
solution  produced  marked  surface  deterioration. 

4.5  Source  material  preparation 

Long,  (ref.  22),  reported  that  GaAs  layers  grown  from  an  undoped  melt 
using  single  crystal  GaAs  source  material  had  lower  impurity  concentra¬ 
tions  than  layers  grown  from  a  similar  melt  using  polycrystalline  source 
material.  He  attributed  the  difference  in  doping  levels  in  the  epilayer 
to  excess  impurities  concentrated  at  the  grain  boundaries  of  the  poly¬ 
crystalline  source,  and  recommended  that  polycrystalline  GaAs  should  only 
be  used  if  the  grain  boundaries  could  be  cut  away  and  discarded.  In 
Melbourne,  polycrystalline  undoped  material  supplied  by  Monsanto  was 
used  as  the  source,  cut  to  fit  within  a  5  mm  x  10  mm  region  in  the  sliding 
boat.  In  Adelaide,  undoped  single  crystal  source  material,  boat  grown 
by  Bell  and  Howell  was  used.  This  was  sawed  into  slices  2  mm  to  3  mm 
thick  and  cut  to  a  size  sufficient  to  cover  the  entire  bottom  of  the 
crucible.  Each  time  a  new  melt  was  prepared,  the  source  block  was 
polished,  cleaned  and  etched  in  a  manner  similar  to  that  used  for  the 
substrate  material,  and  loaded  into  the  furnace  as  rapidly  as  possible. 
Each  melt  mixture  contained  a  large  excess  of  solute,  and  could  thus  be 
used  for  several  successive  growths.  According  to  equation  (l),  a  30  gm 
melt  is  fully  saturated  by  0.42  gm  of  gallium  arsenide  at  700UC  or  by 

1.5  gm  at  800  G.  The  loss  in  weight  of  the  solute  block  after  single  runs 
at  these  temperatures  was,  in  fact,  very  close  to  these  calculated  values. 

4.6  Melt  (solvent)  solution 

As  the  solvent  gallium  was  of  very  high  purity  when  purchased,  no 
further  cleaning  was  attempted.  The  Ga  was  known  to  have  a  low  copper 
content,  which  is  important  since  copper  will  diffuse  rapidly  into 
epitaxial  layers  and  can  cause  undesired  intrinsic  layers. 
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Initially,  gallium  of  seven  9's  purity  from  Alusuisse  Metals  was  used, 
with  30  gm  batches  supplied  in  plastic  bottles.  To  minimise 
contamination  from  the  plastic,  the  bottles  were  refrigerated  until 
immediately  prior  to  use.  When  required,  the  bottles  were  warmed  and 
the  liquid  contents  (melting  point  29.8°C)  poured  over  the  source 
material,  the  contents  of  a  30gm  bottle  covering  the  block  to  a  depth 
of  about  7  mm  in  the  Adelaide  SSM  apparatus.  Because  so  many  different 
melts  were  being  used  for  different  doping  experiments  and  temperature 
profiling,  as  well  as  for  growth  of  undoped  layers,  gallium  of  six  9's 
purity  supplied  by  Cerac  was  decanted  from  a  1000  gm  bottle  into  batches 
of  30  gm  then  stored  and  used  as  detailed  above. 

In  the  TTT  apparatus,  sufficient  liquid  metal  was  squeezed  from  the 
plastic  storage  bottle  to  fill  the  small  bin  and  the  bottle  was  then 
resealed  and  refrigerated  until  the  next  use.  On  each  occasion  only  a 
small  amount  of  solvent  was  used,  of  the  order  of  1  gm  per  growth  run. 

4.7  Impurity  dopants 

The  attempts  to  grow  epilayers  doped  with  a  specific  concentration 
of  zinc  and  silicon  which  were  carried  out  using  the  SSM  were  not  very 
successful,  but  for  completeness  are  reported  below. 

4.7.1  Zinc  doping 

This  group  II  element  acts  as  an  acceptor  and  should  form  p 

1 9 

doped  layers.  For  zinc  concentrations  less  than  about  2  x  10  cm 
essentially  all  the  zinc  atoms  are  ionized  at  room  temperature  and 
so  there  is  very  little  neutral  (uncompensated)  zinc.  The 
liquidus  isotherms  of  Jordan  (reference  24)  were  used  together  with 
the  ternary  phase  diagram  Ga-As-Zn  of  Panish  (reference  25)  to 
calculate  the  mass  of  dopant  Zn  which  should  be  added  to  a  Ga  melt 
in  order  to  achieve  specific  dopant  levels  in  the  epitaxial  layer. 
However,  when  zinc  was  added  to  the  melt,  a  grey- black  deposit 
formed  on  the  quartz  furnace  tube  at  the  operating  temperature. 

This  occurred  for  both  powdered  and  solid  zinc  and  presumably 
results  from  the  high  vapour  pressure  of  zinc  (about  250  torr 
at  800  C).  The  measured  high  resistivity  of  the  epilayer  suggested 
that  little  zinc  had  been  incorporated  into  the  growing  layer,  and, 
indeed,  examination  with  DRCS  X-ray  equipment  could  detect  no  zinc 
at  all. 

Some  interesting  Ga-Zn-As  compounds  with  definite  crystallinity 

were  formed  as  long  silver-grey  threads  elsewhere  in  the  cooler 

parts  of  the  tube.  The  phase  diagram  suggests  that  this  could  be 

either  Zn  As0  or  the  higher  melting  point  Zn„As0.  However,  the 
2  3  2' 

material  has  not  yet  been  analysed. 

4.7.2  Silicon  doping 

Since  the  group  IV  element  Si  is  an  amphoteric  dopant  whose 
behaviour  is  controlled  essentially  by  the  As  vapour  pressure,  it 
can  be  used  for  either  n  or  p  doping.  A  GaAs  crystal  grown  from  a 
melt  close  to  the  stoichiometric  composition,  such  as  our  seed 
crystal,  becomes  n  type  since  Si  is  preferentially  incorporated  on 
Ga  sites  and  acts  as  a  donor.  However,  under  reduced  As  pressure, 
as  occurs  at  lower  temperatures  or  in  a  Ga-rich  melt.  Si  is 
predominantly  incorporated  as  an  acceptor  on  As  sites  and  so  a  p 
type  epilayer  can  be  grown. 
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In  the  experiments  conducted  by  many  different  research  groups,  a 
great  spread  has  been  observed  for  the  temperatures  at  which  the 
conductivity  type  changes  from  n  to  p  for  different  silicon 
concentrations  in  the  melt,  and  correspondingly  large  variations 
in  the  doping  levels  reported  (references  26  -  31 ). 

The  ternary  phase  diagram  Ga-As-Si  proposed  by  Panish 
(reference  32)  together  with  the  liquidus  data  of  Hicks  and  Green 
(reference  28)  was  used  to  predict  the  dopant  concentration  in 
the  epilayer.  In  the  experiments  reported  here,  silicon  of  six  9's 
purity  was  used  without  any  additional  cleaning.  When  0.85  atomic 
percent  of  silicon  was  added  to  a  Ga  melt  saturated  with  GaAs, 
a  smooth  and  shiny  surface  layer  was  grown  at  800°C. 

However,  a  "film"  of  whitish  material  appeared  all  over  the  melt 
surface  and  between  the  melt  and  the  crucible  wall.  This  film  was 
present  at  the  growth  temperature,  and  did  not  appear  only  upon 
cooling  to  room  temperature.  Examination  of  the  melt  showed  that 
less  than  0.08  atomic  percent  Si  had  been  dissolved.  The 
solubility  limit  at  800°C  (reference  33)  should  allow  the 
dissolution  of  up  to  5  atomic  percent  of  silicon,  so  that  all  the 
added  Si  in  the  present  case  should  have  dissolved  and  been 
incorporated  into  the  grown  layer.  A  further  attempt  to  dissolve 
Si  in  Ga  was  made  after  cleaning  the  Si  for  15  seconds  in  a  HF:HU0  = 
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1.3  "white  etch".  Again  very  little  Si  was  dissolved  and  a  similar 
white  film  was  produced.  Although  the  grown  layer  had  a  mirror-like 
surface  finish,  the  junction  was  far  from  planar  and  a  great  deal  of 
further  experimentation  undoubtedly  needs  to  be  carried  out  with  the 
various  dopants. 

4.8  Loading  procedure 

After  preparation  and  final  cleaning,  the  source,  substrate  and  melt 
were  loaded  under  laminar  flow  conditions  and  placed  in  the  furnace  hot 
zone,  ensuring  minimum  exposure  to  laboratory  air  during  the  loading 
operations  by  loading  with  a  N2  atmosphere  flowing  through  the  tube.  The 

tube  was  evacuated  and  pure  hydrogen  admitted.  Once  H2  was  flowing  through 

the  tube  at  atmospheric  pressure,  the  furnace  elements  were  turned  on  and 
the  melt  and  source  allowed  to  equilibriate  for  24  hours.  Several  runs 
were  successfully  carried  out  after  only  14  hours  equilibration,  but 
whenever  a  fresh  melt  was  used,  patchy  and  uneven  growths  were  obtained 
unless  equilibration  times  of  the  order  of  24  hours  were  allowed. 

To  initiate  growth  after  the  equilibration  time,  the  substrate  was 
placed  in  position  in  the  melt,  either  by  sliding  in  TTT,  or  rotating 
and  dipping  in  the  SSM  crucible  arrangement,  and  the  desired  temperature 
cycle  begun.  For  best  results,  one  further  important  step  was  carried 
out  in  the  case  of  the  TTT  experiments.  The  melt  temperature  was  raised 
by  1  to  2  C  for  a  short  time,  to  melt  back  the  substrate  surface  about 
1  /im  and  so  remove  any  chemical  impurities  or  mechanical  damage  remaining 
on  the  surface.  This  ensured  a  clean  interface  on  which  to  grow  the 
epilayer. 

Upon  conclusion  of  the  growth  cycle,  the  substrate  was  removed  from 
the  melt  but  not  unloaded  from  the  H2  atmosphere  until  the  furnace  had 

cooled  to  room  temperature.  Frequently  some  unwanted  Ga  droplets  adhered 
to  the  substrate  surface  and  caused  additional  growth  during  the  cooling 
period.  The  sliding  boat  technique  used  a  scraper  arrangement  to  minimise 
this  adhesion  and  the  consequent  undesired  growth,  but  the  dipping 
technique  relied  only  on  gravity  and  surface  tension  to  remove  excess  melt 
from  the  substrate  surface. 


-  13  - 


ERL-0078-TR 


5.  PROPERTIES  OF  THE  GROWN  EPILAYERS 

For  acceptable  yields  using  epitaxial  growth  in  device  technology,  it  is 
essential  that  the  electrical  and  structural  properties  of  the  grown  layers 
are  well  controlled  and  reproducible.  The  major  properties  of  significance 
for  device  fabrication  are  surface  morphology,  layer  thickness,  junction 
planarity,  doping  level,  and  doping  profile. 

5.1  Surface  Morphology 


Despite  great  care  in  substrate  preparation,  macroscopic  imperfections 
whose  causes  are  not  well  understood,  frequently  occur  in  LPE  layers. 

These  imperfections  can  cause  electrical  inhomogenities  in  injection 
lasers  and  Gunn  oscillators  and  inhibit  their  operation.  The  most  commonly 
reported  surface  morphologies  on  LPE  GaAs  are  surface  terraces,  ripples, 
and  cusps,  whose  shape  and  size  may  be  dependent  on  substrate 
orientation.  Bauser  (ref  34)  claims  that  substrates  must  be  oriented  to 
better  than  5  minutes  of  arc  to  avoid  terrace  structure,  whereas  generally 
available  substrate  material  is  cut  parallel  to  (100)  planes  only  within 
1  or  2  .  The  origin  of  these  surface  imperfections  has  also  been 
attributed  to  substrate  defects,  dislocations,  preparation  residues  and 
scratches,  arsenic  evaporation  etc.,  interface  attachment  (substrate 
orientation,  nuclei  formations),  and  growth  instabilities  (constitutional 
supercooling,  thermal  fluctuations). 

It  had  been  believed  that  constitutional  supercooling  (the  depression 
of  the  liquidus  tempera ture  of  the  melt  at  the  melt-solid  interface 
caused  by  depletion  of  the  solute  during  growth)  could  result  in  an 
unstable  growth  interface  and  cause  spurious  nucleation  of  solid 
throughout  the  melt.  Such  nucleation  in  front  of  the  substrate  surface 
might  result  in  metallic  inclusions,  non-uniform  growth,  and  cellular  growth 
structures  in  the  epitaxial  layer.  However  it  is  now  known  that  there  is 
considerable  difficulty  in  nucleating  the  growth  of  GaAs  from  Ga-rich 
solutions.  Saturated  solutions  can  routinely  be  supercooled  by  greater 
than  5  C  without  any  spontaneous  nucleation  occurring,  and  they  can 
withstand  appreciable  supercooling  even  when  in  contact  with  a  GaAs 
substrate.  Thus  if  a  solution  which  is  just  saturated  (i.e.  at  its 
liquidus  temperature)  is  placed  in  contact  with  a  substrate  at  the  same 
temperature  (in  exact  equilibrium) ,  the  nucleation  will  tend  to  occur  only 
on  a  relatively  small  number  of  favourable  sites.  Because  of  the  high 
mobility  of  molecules  in  solution,  these  "islands"  will  tend  to  grow  at 
the  expense  of  later  forming  sites,  and  gradually  form  a  structured 
surface.  Solomon  (ref  19)  and  Longo  (ref  35)  have  claimed  that  the 
growth  will  remain  discontinuous  until  a  thickness  greater  than  1  ^u  has 

formed.  Supercooling  the  solution  before  inserting  the  substrate  (not 
easily  done  with  the  steady  state  method)  will  increase  the  driving  force 
for  nucleation  and  overcome  the  surface  barrier  to  growth  which  is  evident 
at  equilibrium,  thus  allowing  the  growth  of  smoother  and  more  uniform 
layers.  By  using  this  technique.  Bell  Telephone  Laboratories  (USA)  and 
Standard  Telecommunications  Laboratories  (UK)  have  in  fact  grown  planar 


layers  as  thin  as  0.04  ^u,  in  multiple  layer  growths  (ref.  36,  37). 

For  the  steady  state  growth  method  used  in  Adelaide,  probably  the 
overriding  effect  on  surface  quality  and  layer  thickness  is  the  incomplete 
drainage  of  the  melt  from  the  film  surface  after  removal  from  the  melt. 
Before  it  is  lifted  from  the  melt  the  substrate  holder  is  rotated  as  far 
as  the  crucible  size  permits,  but  this  is  still  substantially  less  than 
the  90  required  for  gravity  and  surface  tension  to  achieve  their  maximum 
drainage  effect. 
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The  surfaces  of  these  layers,  with  the  exception  of  "tide-marks"  sometimes 
noticed,  are  generally  smooth  and  highly  reflective  to  the  naked  eye,  but 
exhibit  pronounced  surface  structure  when  viewed  under  Nomarski 
interference  contrast  on  an  optical  microscope.  (Figures  2(a)  and  2(b) 
for  two  separate  substrates  grown  by  SSM  show  how  great  the  variation  in 
surface  morphology  is  over  each  individual  substrate). 

5.2  Doping  level  and  doping  profile 

Once  the  thickness  of  the  grown  layer  has  been  determined  (using  the 
techniques  described  in  Section  5.4)  the  average  free  carrier 
concentration  over  that  portion  of  the  grown  layer  under  examination 
can  be  determined  by  measurements  of  the  Hall  effect.  Measurements  of 
the  junction  capacitance  can  also  be  used  to  predict  doping  levels,  and 
can  in  addition  indicate  whether  the  grown  layer  forms  an  abrupt  junction 
or  whether  diffusion  of  carriers  has  taken  place,  which  would  result 
in  a  graded  junction.  Although  these  techniques  were  used  at  the 
Telecom  Laboratories  on  some  gallium  arsenide  samples  having  p-n  junctions 
formed  by  diffusion,  no  such  analyses  for  epitaxial  samples  were  at  that 
time  attempted. 


5.3  Junction  Planarity 


If  the  substrate  was  positioned  directly  above  the  melt  during  the 
equilibration  period  where  it  was  exposed  to  volatile  material  driven  off 
from  the  melt,  the  layers  often  grew  unevenly.  Patchy  growths  of  this 
nature  are  seen  in  figures  3(a)  and  3(b),  where  a  series  of  mesas  have 
been  formed  on  the  substrate.  This  type  of  behaviour  was  most  pronounced 
whenever  a  new  melt  was  used  for  the  first  time. 

When  the  substrate  was  protected  from  the  vapours  until  immediately 
prior  to  the  commencement  of  growth,  good  flat  junctions  over  most  of  the 
substrate  were  routinely  achieved.  The  junctions  were  made  visible  by 
cleaving  the  substrate  along  ( 1 1 0 )  planes  and  stainihg  the  cleaved 
section  with  a  solution  of  1HF:3HN0^:4H20  for  15-40  seconds.  Differences 

in  electron  densities  in  the  differently  doped  regions  affect  the  rate  at 
which  GaAs  reacts  with  the  staining  solution,  and  the  difference  in 
etching  rates  causes  the  junction  to  become  visible  under  a  microscope. 
(This  method  requires  different  doping  levels,  but’ not  necessarily 
different  carrier  types  -  although  p-n  junctions  are  better  delineated 
than,  for  example,  undoped  -n  or  n-n+  junctions). 

The  typically  flat  junction  formed  on  a  substrate  is  shown  in  figure  4. 
However  edge  effects  of  some  sort  are  present,  as  the  junctions  are 
frequently  not  flat  near  the  edges  of  the  substrate,  and  odd  shaped 
overgrowths  are  also  sometimes  observed  at  the  corners.  The  very  irregular 
junction  obtained  after  the  attempted  growth  of  a  silicon  doped  epilayer 
(described  in  Section  4.7.2)  is  shown  in  figure  5. 

5.4  Epilayer  thickness 


The  thickness  of  a  grown  layer  was  measured  after  the  cleaving  and 
staining  procedure  mentioned  in  Section  5.3«  To  view  the  delineated 
junctions  and  estimate  the  thickness  of  the  grown  layers,  a  Zeiss  KK08 
reflected  light  microscope  was  used  with  a  Nomarksi  differential 
interference  contrast  attachment  and  a  magnification  of  128X.  Calibration 
of  the  eyepiece  scale  yielded  a  value  of  6.8  ^u  per  small  division  at  this 

magnification,  so  that  the  layer  height  could  be  resolved  to  within  a  few 
microns  at  each  position.  Figure  6  is  a  schematic  representation  of  the 
thickness  determined  in  this  manner  for  one  sample  where  steady  state 
epitaxial  growth  was  undertaken  over  a  period  of  2  hours  at  nominally 
800°C  (the  furnace  tog  element  was  maintained  at  770  C  while  the  bottom 
was  maintained  at  820  C). 
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It  is  evident  that  there  is  a  Great  variation  in  the  height  of  the  epilayer 
over  this  substrate,  and  similar  variations  were  found  in  most  samples 
grown  at  Adelaide  using  the  steady  state  method.  However,  as  is  more  fully 
reported  by  the  present  author  in  reference  16,  the  precise  variations  in 
height  of  the  grown  epilayers  differ  for  each  sample.  Long  (reference  15) 
in  the  only  published  work  using  the  SSM,  has  reported  that  2:1  variations 
in  layer  height  across  the  surface  of  a  single  substrate  are  common,  and 
his  observations  are  consistent  with  the  present  results. 


6.  DISCUSSION 


The  variable  epilayer  thickness  over  a  given  substrate,  and  the  non¬ 
reproducibility  of  results  from  sample  to  sample,  strongly  suggest  that 
processes  other  than  simple  diffusion  of  gallium  arsenide  through  the  melt  from 
source  to  substrate  are  taking  place.  The  layers  grown  by  the  Steady  State 
method  are  in  general  much  thicker  than  predicted  by  the  diffusion  equation 
developed  in  Section  3,  which  suggests  that  the  transfer  of  solute  material 
to  the  growth  interface  is  greater  than  for  diffusion-limited  "steady  state" 
conditions.  In  addition,  under  Nomarski  differential  interference  contrast 
on  an  optical  microscope  the  grown  layers  exhibit  flow  marks  and  irregular 
wave-like  patterns  on  their  surfaces.  The  irregular  surface  morphology  might 
be  explicable  in  terms  of  unstable  conditions  created  at  the  growth  interface 
by  fluid  motion.  Any  point  on  the  surface  of  the  growing  layer  would  then  be 
subjected  to  a  supersaturation  which  varies  with  time. 

If  fluid  motion  is  occurring  in  the  melt,  then  kinetic  processes  such  as 
convection  will  be  operative.  Convective  forces  could  result  from  horizontal 
temperature  variations  within  the  melt,  which  will  develop  if  the  furnace 
temperatures  are  not  strictly  uniform  in  the  longitudinal  direction.  In  fact, 
when  the  Marshall  split  furnace  used  for  SSM  was  carefully  profiled 
(reference  16),  the  temperature  of  the  top  element  at  a  horizontal  distance 
2  cm  from  the  mid-gosition  was  found  to  be  15°C  greater  than  the  set 
temperature  of  770  C  at  the  mid-position,  while  the  bottom  element  was 
reasonably  uniform  over  this  region.  These  large  horizontal  variations  in 
the  temperature  differential  AT  across  the  furnace  elements  occur  in  the 
region  where  the  substrate  is  positioned  during  a  growth  run,  and  thus  would 
be  expected  to  produce  a  horizontal  variation  in  temperature  within  the  melt 
itself.  Such  variations  were  indeed  observed  when  a  differential  thermocouple 
probe  was  positioned  actually  within  the  melt  (reference  16). 

Because  the  solubility  of  gallium  arsenide  in  gallium  is  temperature 
dependant,  then  the  existence  of  horizontal  as  well  as  vertical  temperature 
gradients  within  the  melt  will  cause  variations  in  solute  concentration  in 
both  horizontal  and  vertical  planes.  Cellular  convection  in  the  melt  could 
be  caused  by  direct  thermal  gradients,  or  else  by  the  solutal  gradients 
formed  by  these  temperature  variations.  In  Appendix  II  we  develop  a  treatment 
in  terms  of  the  Raleigh  criteria  (references  15,  38)  to  determine  the 
conditions  for  the  onset  of  convection  produced  by  either  temperature  gradients 
Ay  or  concentration  gradients  AC  .  The  results  in  Tables  II. 1  and  II. 2  for 
Ax  Ax 

the  experimental  configuration  used  in  the  SSM  show  that  the  direct 
contribution  of  the  temperature  gradient  to  solution  instability  is  much  less 
than  the  contribution  of  the  corresponding  temperature-induced  concentration 
gradient.  Thus  solute  gradients  rather  than  thermal  gradients  are  the  driving 
force  which  cause  convection  to  develop  within  the  melt. 

The  Raleigh  criteria  predict  that  solution  heights  less  than  a  critical 
value  are  necessary  to  ensure  that  convection  currents  do  not  develop. 
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If  we  label  the  critical  values  Ax^,  and  Ax^  to  represent  temperature-induced 

and  solute-induced  effects  respectively,  then  the  experimental  conditions  at 
800  C  with  /  AT  \  =  4.6°C/cm 

V  Ax  / 

yield 

Ax^  =  10.9  mm 

Axq  =  1 .3  mm 

In  the  majority  of  growths  carried  out  using  the  SSM,  a  30  g  gallium  melt 
was  used.  This  quantity  of  gallium  provides  a  solution  depth  of  7  mm,  while 
the  substrate  is  positioned  at  heights  varying  from  2  mm  to  5  mm.  Thus  it  is 
clear  that  convection  will  be  produced  in  this  configuration,  resulting  in 
the  non-uniform  grown  layers  which  were  in  fact  observed. 

Since  this  convection  provides  a  mechanism  for  non-uniform  growth,  it  must 
be  eliminated,  or  at  least  controlled,  to  obtain  uniform  thickness  and  smooth 
surfaces  for  the  grown  epitaxial  layers.  The  calculations  in  Appendix  II 
demonstrate  that  very  small  heights  (~  1  mm)  are  necessary  to  inhibit 
convection  effects  when  the  GaAs  substrate  is  placed  towards  the  top  of  the 
solution.  However  if  the  substrate  is  placed  below  the  saturated  melt,  and 
the  temperature  gradient  is  reversed  (so  that  the  substrate  is  still  at  a 
lower  temperature  than  the  source,  and  the  solution  is  everywhere  warmer  than 
the  substrate)  then  the  depletion  of  solute  which  occurs  adjacent  to  the 
substrate  will  result  in  a  stabilising  density  gradient  (ref  11,  15).  There 
will  be  little  tendency  for  solute-induced  convection  to  develop,  and 
temperature-induced  convection  should  not  develop  until  A  x  >  io  mm.  Using  the 
SSM  with  the  substrate  below  the  melt,  the  Cornell  group  (reference  15)  report 
that  uniform  epitaxial  layers  with  heights  close  to  those  predicted  by  the 
diffusion  -  limited  conditions  of  equation  (4)  were  observed. 

However  for  the  type  of  devices  of  interest  to  DECS,  the  temperature 
transient  technique  using  a  sliding  boat  arrangement  with  small  melt  size 
and  with  the  substrate  positioned  under  the  melt  bins  is  more  attractive  than 
the  steady  state  method  using  large  melt  volumes  in  a  relatively  large  crucible. 
The  sliding  boat  technique  permits  the  growth  of  several  consecutive  epitaxial 
layers  merely  by  sliding  the  substrate  from  under  one  bin  to  under  another.  The 
layers  may  be  of  different  compositions  depending  on  the  composition  of  each 
melt,  and  the  thickness  of  each  layer  is  controlled  by  the  time  the  substrate 
is  in  contact  with  the  corresponding  melt.  Carry-over  of  solution  from  bin  to 
bin  may  be  minimised  by  suitable  construction  of  the  slider. 

On  the  other  hand,  using  the  steady  state  method  where  the  substrate  is 
positioned  towards  the  top  of  the  melt  solution  causes  convection  effects  in 
the  melt  and  consequent  non-uniformities  to  develop  in  the  growing  layer.  When 
the  substrate  is  mechanically  moved  from  one  crucible  to  another  in  the  attempt 
to  grow  multiple  layers,  there  is  likely  to  be  a  greater  carry-over  of  solution 
than  in  the  sliding  boat  technique.  More  importantly,  since  the  surface  of  the 
first  epilayer  is  not  flat,  subsequent  epitaxial  layers  will  not  form  flat 
junctions  and  the  non-uniform  multiple  layers  will  not  be  of  use  for  making 
practical  devices.  Although  more  uniform  epitaxial  layers  will  be  grown  if 
the  substrate  is  held  below  the  saturated  melt  rather  than  some  distance  above 
the  source  material,  such  a  configuration  is  difficult  to  achieve  with  the 
present  apparatus.  In  addition  there  is  then  no  suitable  way  of  moving  the 
substrate  to  another  crucible,  so  that  multiple  layer  growths  are  not  possible. 
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7.  CONCLUSION 


Epitaxial  growth  of  gallium  arsenide  from  the  liquid  phase  is  discussed 
for  both  the  temperature  transient  technique  and  the  steady  state  method,  and 
the  properties  of  epitaxial  layers  grown  under  a  variety  of  experimental 
conditions  and  in  two  different  laboratories  are  presented.  The  observed 
heights  of  the  grown  layers  are  compared  with  mathematical  models  which  were 
developed  for  Doth  growth  processes.  Although  reasonable  agreement  is  found 
between  calculated  and  observed  layer  heights  for  the  temperature  transient 
technique,  a  very  large  discrepancy  is  found  when  the  present  form  of  the 
Steady  State  method  is  employed.  The  grown  layers  are  much  thicker  than  would 
be  expected  if  the  growth  was  limited  by  diffusion  of  solute  through  the  melt, 
in  addition  to  being  of  non-uniform  thickness.  The  non-uniformity  suggests 
that  convection  effects  might  be  present  in  the  growth  solution,  a  possibility 
which  is  discussed  at  some  length.  It  is  demonstrated  that  for  the  apparatus 
used  here  for  the  steady  state  method,  the  temperature  non-uniformities  which 
are  produced  within  the  growth  solution  are  indeed  sufficient  to  allow 
convective  forces  to  develop.  The  relative  advantages  of  the  two  growth  methods 
are  discussed,  and  it  is  concluded  that  the  temperature  transient  technique 
using  a  multiple-bin  sliding  boat  is  more  suitable  for  forming  the  type  of 
devices  of  interest  to  DRCS.  The  experience  gained  using  the  different 
techniques,  and  the  realization  of  the  critical  importance  of  temperature 
uniformity  and  stability,  have  led  to  the  design  of  an  improved  system  for 
liquid  phase  epitaxial  growth  at  DRCS.  The  very  satisfactory  initial  results 
which  have  been  obtained  with  this  new  system  encourage  the  belief  that  it  can 
fulfil  the  requirements  of  an  "ideal"  epitaxial  system  and  enable  the  growth 
of  uniform  and  controllable  layers  suitable  for  the  manufacture  of  practical 
electro-optic  devices. 
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APPENDIX  I 

1.1  Slope  oi‘  the  Liquidus  Curve 

The  slope  of  the  liquidus  curve  is  obtained  from  Hall's  solubility  data 
(reference  14)  which  receives  independent  verification  from  Hsieh 
(reference  5)  in  a  limited  temperature  range  near  800  C. 


TABLE  1.1 

t(°c) 

XAS  (T) 

xL (T) 

dX  (t)  (mole  fraction 

As  0 

(mole  fraction/  C) 

(mole  fraction/  C) 

dT  per  C  ) 

700 

800 

900 

6.50  x  10“p 

2.32  x  10  p 

5.69  x  10 

6.59  x  10"p 

2.43  x  10 

6.75  x  10 

0.90  x  10“J 

2.78  x  10“T 

7.98  x  10  4 

1.2  The  Diffusion  Coefficient  D 

-4  -5 

Diffusion  coefficients  for  most  liquids  lie  in  the  range  10  to  10 

2 

cm  /sec,  with  some  temperature  dependence  observed.  No  direct  experimental 
measurements  of  the  diffusion  of  As  or  GaAs  in  a  Ga  solvent  have  been 
published.  Instead,  estimates  for  D  have  been  made  by  comparing 
experimental  growth  rates  with  the  predictions  of  various  mathematical 
theories  (references  5,6,7,39,40),  and  obtaining  "best  fit"  values  for  D. 
There  is  a  large  variation  in  the  values  so  obtained,  but  based  on  the 
published  results,  we  consider  the  values  of  Table  1.2  appropriate  to  the 
present  experimental  conditions. 


TABLE  1.2 


T  (°C) 

2 

D  (cm  /sec) 

700 

800 

900 

3  x  10"^ 

5  x  10  ? 

10  x  10  5 

1.3  Solution  Temperature  Gradient 

Temperature  probing  within  the  melt  was  undertaken  and  is  described  in 
detail  in  reference  16. 

To  make  these  measurements,  a  thermocouple  made  from  0.025"  diameter 
chromel  and  alumel  wires  was  encased  in  thin  quartz  and  positioned  with 
its  "hot"  junction  at  the  bottom  of  a  probe  which  had  dimensions  similar 
to  the  substrate  holder  used  in  SSM  growths.  The  probe  was  dipped  into 
the  melt  and  then  raised  or  lowered  a  certain  distance  Ax  as  determined 
by  a  simple  screw  arrangement  external  to  the  tube.  This  was  very  crude, 
but  allowed  an  order-of-magnitude  determination  of  the  temperature  gradient. 
The  repeatability  of  measurements  from  day  to  day,  or  even  in  successive 
cycles  of  the  probe,  was  poor.  However,  the  results  indicated  that  the 
solution  temperature  gradient  was  almost  independent  of  the  actual  furnace 
temperature  and  largely  determined  by  the  temperature  differential  between 
the  top  and  bottom  elements.  For  an  element-to-element  differential 
iT  =  50°C,  the  solution  temperature  gradient  was  found  to  be  of  the 
order  of 

4.2°C/ cm  at  700°C  and 
4.6°C/cm  at  800°C. 


-  23  - 


ERL-0078-TR 


APPENDIX  II 

Convection  Effects  in  Fluids 


In  a  two  component  system,  such  as  a  gallium  melt  containing  dissolved 
arsenic,  there  are  two  major  processes  which  could  contribute  to  natural 
convection.  Instabilities  may  be  created  by  either  thermal  gradients  or 
concentration  gradients  in  the  solution.  These  potential  instabilities  are 
opposed  by  the  viscosity  of  the  fluid,  which  provides  a  damping  force. 
However  instability  and  consequent  cellular  convection  will  occur  whenever 
a  critical  temperature  gradient  or  a  critical  solute  gradient  (which  may  be 
temperature  induced)  is  exceeded.  These  critical  values  depend  on  the 
particular  experimental  situation  and  geometry. 

Raleigh  (ref  38)  originally  showed  that  a  fluid  subjected  to  a  purely 
vertical  temperature  difference  AT  across  a  depth  of  fluid  Ax  can  be 
characterised  by  a  Raleigh  number  R^,  where 


Rt  =  g  AT  y 


( An3 

K  i, 


K  T) 


(Ax)4 


(II.  1  ) 


and  V  =  coefficient  of  volumetric  expansion  (=  9.7  x  10  '’/°C  at  900°C  for  Ga) 

2 

K  =  thermal  diffusivity  (=  0.3  cm  /s) 

V  =  kinematic  viscosity  (=  1.2  x  10  ^  Stokes  at  800°C) 

2 

g  =  980  cm  /s 


(The  values  of  the  parameters  quoted  above  were  obtained  from  reference  15). 

Instability  with  respect  to  the  thermal  gradient  occurs  when  the  Raleigh 

number  R„  exceeds  a  certain  critical  value  R  .  For  a  fluid  bounded  by  two 
i.  c 

solid  surfaces,  this  critical  value  is  Rc  =  (1700  +_  140).  Thus  instability 
will  set  in  when 

Rt  >(1700  +  140). 

If  the  solute  is  less  dense  than  the  solvent,  it  is  possible  to  get  density 
inversion  within  the  melt  when  the  substrate  is  at  a  lower  temperature  than 
the  source  and  is  located  above  the  source  (ref  11 ).  An  expression  similar  to 
equation  ( II. 1 )  may  be  derived  (ref  11,  15,  41 )  for  differences  in  solute 
concentration  AC  over  a  distance  Ax  in  the  fluid. 


where  /3  =  volumetric  change  for  unit  concentration  change 

D  =  mass  diffusivity  (=  5  x  10  ^  cm^/s  at  800°c). 

Instability  with  respect  to  the  density  gradient  produced  by  a  single  solute 
sets  in  when 


Rg  >1770. 
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When  temperature  gradients  and  solute  gradients  are  both  acting,  as  in  the 
present  experimental  situation,  then  the  individual  effects  will  be  superimposed. 
Because  the  thermal  diffusivity  K  is  much  greater  than  the  mass  diffusivity  D, 
the  critical  Raleigh  number  will  be  exceeded  more  readily  as  a  result  of 
concentration  variations  than  as  a  result  of  temperature  differences.  When 
the  critical  value  is  exceeded  by  even  a  small  value,  steady  convection  cells 
are  set  up.  If  the  critical  value  is  exceeded  by  an  order  of  magnitude  or 
more,  then  the  convection  becomes  turbulent,  and  irregular.  In  fact,  a 
solutally  convecting  solution  is  likely  to  develop  a  changing  wave  pattern 
on  the  interface,  which  for  thick  layers  may  develop  into  a  cellular  pattern 
reminiscent  of  the  effects  of  constitutional  super-cooling  in  metals  (ref  11 ). 

Using  the  critical  Raleigh  numbers,  it  is  possible  to  calculate  as  in 
section  II.  1  and  II. 2  below  the  approximate  experimental  conditions  under 
which  stable  growth  should  result. 


II.  1  Temperature-induced  instabilities 


The  stability  criterion 


£L 

K7 


(Ax)4 


<  1700 


predicts  that  solutions  of  depth  Ax  will  be  stable  against  convection 
for 


(  Ax)4 


<  6.45 


at  800°C 


Thus  stability  with  respect  to  thermal  gradients  will  exist  for  values 
of  Ax  less  than  the  values  in  Table  II. 1. 


TABLE  II.  1 


(  at)  o_  / 
l  ~a —  y  c/cm 
\Ax/ 

Ax  cm 

0.25 

2.25 

0.5 

1  .89 

1.0 

1  .59 

2.0 

1  .54 

4.0 

1.15 

5.0 

1 .07 

10.0 

0.90 

/  A.T  \  o 

In  the  present  experimental  situation  where  1— -  )=  4.6  C/cm,  it  is  necessary 
to  have  Ax  < 1 .09  cm  to  prevent  convection  developing  in  the  solution. 


II. 2  Solute-induced  instabilities 


The  stability  criterion  for  solute  concentration  gradients 


is  more  difficult  to  analyse  than  that  for  temperature  gradients.  Wo 
experimentally  determined  values  exist  for  jS ,  the  change  in  volume  of  a 
solution  for  unit  change  in  solute  concentration. 
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A  complicated  formulation  by  Tiller  (ref.  41 )  may  be  simplified  (ref.  15) 
for  the  Ga-As  binary  system  to  the  form 


+ 


where  is  the  concentration  of  As  at  the  solid-liquid  growth  interface. 

The  density  of  free  As,  p  (measured  at  20°C)  is  greater  than  the 

density  p ^  of  the  solvent,  which  would  tend  to  discourage  solute-induced 

instabilities.  However,  there  are  reasons  to  believe  that  As  does  not  exist 
freely  in  the  solution,  but  that  some  chemical  bonding  Ga-As  occurs 
(ref.  42-44).  If  the  solute  As  does  exist  as  a  Ga-As  complex,  then  we  should 
replace  p ^  by  and  by  MGaAg,  in  calculating  a  value  for  /?  . 

Using  this  assumption,  a  value  for  j3  may  be  estimated  as  2.1  x  10~^/(g  cm-'’). 

The  stability  criterion  then  becomes 


i.e. 


(  Ax)4  <  15  x  10  4 


This  predicts  that  stability  with  respect  to  solute  concentration  gradients 
will  exist  only  for  values  of  Ax  less  than  the  values  in  Table  II. 2. 


TABLE  II. 2 


necessary  to  have  A x <  0.13  cm  to  prevent  convection  developing. 
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Figure  1 


GaAs 


(a)  Sliding  boat  used  for  the  Temperature  Transient  Technique 


Pyrolytic  graphite 

substrate  holder  Quartz  control  rod 


crucible 


(b)  Crucible  arrangement  used  for  the  Steady  State  Method 

Figure  1 .  Schematic  of  apparatus  used  for  liquid  phase  epitaxial 
growth  of  a  single  layer 
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Figure  2(a) 
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Figure  2(b) 


Scale :  71  \i/  cm 

Figure  2(b).  Various  surface  morphologies  observed  on  Sample  9 
(Magnification  140X; 
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Figures  4&5 


Junction 


Figure  4.  Highly  uniform  junction  formed  on  a  GaAs  substrate 


Si  doped  epilayer 
Irregular  junction 


Figure  5.  Irregular  junction  formed  after  attempts 
to  dope  the  epitaxial  layer  with  silicon 
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Figure  6 


Identification  mark 


Figure  6.  Topical  example  of  the  variation  in  thickness  of  the 
epitaxial  layers  grown  using  the  Steady  State  Method. 
The  sample  shown  here  was  grown  over  a  period  of 
2  hours,  at  a  nominal  temperature  of  800°C  and  with 
a  temperature  differential  of  50  C  applied  across 
the  furnace. 
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